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HIGHLIGHTS 


•  The  thermal  conductivity  of  the  membrane  was  comparable  in  both  directions. 

•  The  through-plane  thermal  conductivity  of  the  membrane  decreases  with  temperature. 

•  The  through-plane  thermal  conductivity  of  the  GDL  decreases  with  temperature. 

•  The  through-plane  thermal  conductivity  of  the  GDL  increases  with  compression. 

•  The  thermal  conductivity  of  the  catalyst  was  comparable  in  both  directions. 
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The  thermal  conductivity  of  the  components  of  the  membrane  electrode  assembly  (MEA)  and  GDL  must 
be  accurately  estimated  in  order  to  better  understand  the  heat  transfer  processes  in  the  proton  exchange 
membrane  (PEM)  fuel  cells.  In  this  study,  an  experimental  investigation  has  been  performed  to  measure 
the  through-plane  thermal  conductivity  and  the  contact  resistance  for  a  number  of  gas  diffusion  layer 
(GDL)  materials.  The  sensitivity  of  these  quantities  to  the  temperature,  PTFE  content  and  micro  porous 
layer  (MPL)  coating  has  been  undertaken.  In  addition,  the  through-plane  thermal  conductivity  of  the 
membrane  has  been  measured  and  reported  as  a  function  of  temperature  and  water  content.  Further,  the 
through-plane  thermal  conductivity  of  the  catalyst  layer  has  been  determined  as  a  function  of  temper¬ 
ature  and  platinum  loading.  It  has  been  found  that  the  through-plane  thermal  conductivity  of  the 
components  of  the  MEA  decreases  when  the  temperature  increases,  and  the  through-plane  thermal 
conductivity  of  the  GDL  is  significantly  lower  than  its  in-plane  thermal  conductivity. 

©  2014  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

The  thermal  conductivity  of  the  components  of  the  membrane 
electrode  assembly  (MEA)  and  GDL  must  be  accurately  estimated  in 
both  directions,  namely  in-plane  and  through-plane  directions,  in 
order  to  better  understand  the  heat  transfer  processes  in  proton 
exchange  membrane  (PEM)  fuel  cells  [1—3].  Many  researchers  use  a 
steady  state  method  to  measure  the  thermal  conductivity  of  gas 
diffusion  layers  (GDLs)  in  the  through-plane  directions  [4,5].  In 
particular,  Vie  et  al.  [6]  were  the  first  research  group  who  attempted 
to  measure  the  thermal  conductivity  of  the  fuel  cell  components.  In 
their  paper,  many  thermocouples  have  been  inserted  in  the  fuel  cell 
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and  the  temperature  gradient  is  measured  at  different  locations.  The 
thermocouples  were  inserted  between  the  gas  diffusion  layers,  the 
catalyst  layers  and  the  membrane.  The  thermal  conductivity  of  the 
E-Tek  ELAT  GDL  and  the  catalyst  layer  was  about  0.2  ±  0.1  W  m-1  K-1. 
However,  these  measurements  were  not  accurate  due  to  the  high 
uncertainty  on  the  locations  of  the  thermocouples  and  due  to  the 
fact  that  the  thermocouples  blocked  some  of  the  active  area  of  the 
fuel  cell.  Khandelwal  and  Mench  [7]  measured  the  through-plane 
thermal  conductivity  of  Toray  and  SIGRACET®  carbon  papers  under 
20  bar  compression  pressure.  Ramousse  et  al.  [8]  estimated  the 
thermal  conductivity  of  a  typical  GDL  to  be  lower  than  the  thermal 
conductivity  of  pure  carbon  samples  but  they  did  not  take  into  ac¬ 
count  the  effect  of  the  compaction  pressure  on  the  thickness  of  the 
sample.  Nitta  et  al.  [9]  measured  the  thermal  conductivity  of  SGL 
10BA  GDL  and  the  thermal  contact  resistance  between  the  GDL  and 
the  graphite  rods.  It  was  found  that  the  values  of  the  GDL  thermal 
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conductivity  obtained  were  almost  4  times  larger  than  those  found 
in  the  literature  and  it  depends  on  the  compression  pressure  on  the 
sample.  They  reported  the  measured  through-plane  thermal  con¬ 
ductivity  of  the  GDL  to  be  about  1.8  ±  0.11  W  itT1  K"1.  Karimi  et  al. 
[10]  determined  the  through-plane  thermal  conductivity  of  Spec- 
traCarb  GDL  experimentally.  The  contact  resistance  between  the 
GDL  and  the  aluminium  apparatus  surface  was  studied  as  a  function 
of  compression  and  PTFE  content.  The  results  showed  that  as  the 
compression  load  increases,  the  thermal  conductivity  of  the  GDL 
increases.  The  obtained  through-plane  thermal  conductivity  of  the 
GDL  was  0.36  and  0.7  W  itT1  Rounder  a  compression  load  0.7  and 
13.8  bar,  respectively.  Burheim  et  al.  11  ]  reported  that  the  through- 
plane  thermal  conductivity  of  dry  GDLs  and  GDLs  that  contain  water. 
The  liquid  water  was  sucked  for  about  30  s  and  the  volume  fraction  of 
the  water  was  calculated  by  comparing  the  weight  of  the  dry  and  wet 
GDLs.  Their  results  showed  that  the  through-plane  thermal  con¬ 
ductivity  of  the  GDL  increases  with  the  presence  of  liquid  water  and 
it  is  0.15  and  1.6  W  m”1  K-1  for  the  dry  and  wet  GDL,  respectively. 

In  this  paper,  an  experimental  setup,  based  on  the  steady-state 
method,  is  developed  to  measure  the  through-plane  thermal  con¬ 
ductivity  of  the  components  in  the  membrane  electrode  assembly 
(MEA)  and  GDL  at  different  operating  temperatures.  The  thermal 
conductivities  of  the  GDLs  are  investigated  as  a  function  of  the  PTFE 
loading,  temperature  and  compression  pressure.  In  addition,  for  the 
present  study  to  be  comprehensive,  the  through-plane  thermal 
conductivities  of  Nation8  membranes  and  catalyst  layers  are 
measured  and  reported  as  a  function  of  the  temperature. 

2.  Materials  and  procedures 

2.2.  Test  apparatus 

An  experimental  apparatus  has  been  developed  to  measure  the 
thermal  conductivity  of  the  various  components  of  the  MEA  and  GDL 
under  steady  state  conditions.  Therefore,  the  formula  employed  to 
estimate  the  thermal  conductivity  is  the  Fourier  law  [12]: 


Ls 

where  As  is  the  cross-sectional  area  of  the  sample,  Ls  is  the  length  of 
the  sample,  ks  is  the  thermal  conductivity  of  the  sample,  and  AT  is 
the  temperature  drop  across  the  sample. 

The  test  apparatus  is  shown  in  Fig.  1.  It  consists  of,  from  top  to 
bottom,  (i)  a  dial  gauge  indicator  to  measure  the  reduction  in  the 
thickness  of  the  sample  under  compression,  (ii)  a  load  cell  which 
records  the  compression  pressure  on  the  sample,  (iii)  the  upper 
steel  flux  metre,  which  contains  3  thermocouples,  (iv)  the  tested 
sample,  and  (v)  the  lower  steel  flux  meter  which  also  contains  3 
thermocouples  whose  temperature  gradient  is  maintained  low  and 
constant  using  a  cooling  system,  see  Section  2.3  for  more  details. 

2.2.  Materials 

The  through-plane  thermal  conductivity  of  the  GDL  is  deter¬ 
mined  for  five  different  SGL  samples  (10AA,  10BA,  10CA,  10DA, 
10EA)  whose  PTFE  loading  are  0,  5, 10,  20  and  30%,  respectively.  In 
addition,  the  effect  of  the  micro  porous  layer  (MPL)  coating  on  the 
through-plane  thermal  conductivity  has  been  investigated  for  two 
different  samples,  namely  10BC  and  10BE,  and  the  thermal  con¬ 
ductivity  of  them  has  been  compared  with  that  of  their  base  sub¬ 
strate,  namely  10BA.  These  GDL  samples,  which  were  provided  by 
the  SGL  Technologies  Gm  bH,  Germany,  are  listed  in  Table  1. 

The  through-plane  thermal  conductivity  of  a  115  Nation8  115 
membrane  (Du  Pont,  USA),  which  is  about  127  pm  thick,  is  also 


Fig.  1.  Configuration  of  the  experimental  set-up  to  measure  the  through-plane  thermal 
conductivity  of  the  MEA  components. 

measured  and  reported.  Furthermore,  the  through-plane  thermal 
conductivity  of  the  catalyst  is  evaluated  with  three  different  plat¬ 
inum  (Pt)  loadings,  namely  0.2,  0.4  and  0.6  mg  cm-2  in  order  to 
investigate  the  effect  of  this  loading  on  the  through-plane  thermal 
conductivity  of  the  catalyst  layer. 

2.3.  Experimental  conditions 

All  measurements  were  performed  under  vacuum  conditions  in 
order  to  eliminate  the  heat  transfer  by  convection.  Moreover,  the 
fixtures  and  the  samples  were  well  insulated  by  using  Rockwool 
insulation  to  minimise  the  heat  loss  in  the  radial  direction  and 
mitigate  heat  transfer  by  radiation.  The  effect  of  the  temperature  on 
the  through-plane  thermal  conductivity  of  all  the  components  in 
the  MEA  and  GDL  was  investigated  in  the  temperature  range 
35-80  °C,  which  is  the  most  likely  operating  temperature  range  of 
PEM  fuel  cells  13].  In  addition,  the  effect  of  the  compression 
pressure  was  investigated  for  the  compression  range  1-20  bar,  in 
which  the  normally-used  compressive  pressure  on  PEM  fuel  cells 


Table  1 

Manufacturers'  specifications  for  the  tested  GDLs. 


GDL 

Thickness  (pm) 

PTFE  loading  (wt.  %) 

10AA 

390 

0 

10BA 

400 

5 

10CA 

400 

10 

10DA 

400 

20 

10EA 

374 

30 

10BC 

415 

23a 

10BE 

367 

50a 

a  PTFE  loading  in  the  MPL  (wt.  %). 
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lies  [14].  The  lower  steel  cylinder  was  maintained  cold  by  using 
EXT-440  Koolance's  and  CHC-122  cooling  block  through  which  the 
cooling  fluid  flows.  All  the  thermocouple  readings,  the  current,  the 
voltage  in  the  circuit  and  the  pressure  readings,  were  controlled 
and  monitored  by  a  LabVIEW  application,  see  Fig.  2. 

It  is  important  to  note  that  every  measurement  in  this  study  was 
repeated  at  least  twice  in  order  to  cheque  for  repeatability.  If  the 
results  were  not  repeatable  then  further  measurements  were  made 
until  repeatable  results  were  obtained. 

2.4.  Methodology 


The  axial  heat  flux  through  the  tested  samples  was  estimated  by 
averaging  the  heat  fluxes  through  upper  and  lower  flux  meters  by 
using  the  following  equations: 


,  h-h 

4up  —  km  T  T 

L\  ~  l3 

(2) 

j  T4  —  Tq 

4 low  —  km  .  , 

l4  ~  l6 

(3) 

4up  +  4low 

qs~  2 

(4) 

where  qs  is  the  heat  flux  through  the  sample,  qup  is  the  heat  flux 
through  the  upper  flux-metre,  qiow  is  the  heat  flux  through  the  lower 
flux-meter,  km  is  the  thermal  conductivity  of  the  standard  material, 
and  Ti,  T2,  T3,  T4,  T5,  and  Tq  are  the  temperatures  measured  by  the  first, 
second,  third,  fourth,  fifth  and  the  sixth  thermocouples,  respectively. 

The  temperature  drop  AT  across  the  test  sample  is  obtained  by 
'joining'  the  temperature  gradients  of  the  lower  and  upper  fixtures, 
as  shown  in  Fig.  3. 


The  total  resistance,  ftrotah  to  the  heat  flow,  which  is  the  sum  of 
‘bulk’  thermal  resistance  and  the  contact  resistance  [12  ,  can  be 
expressed  using  the  Fourier  equation  as  follows: 

^Total  =  ~7T~  =  2^/-/- sample  3“  ^th  (5) 

qs 

where  Rt h  is  the  sample  thermal  resistance,  i?H-sampie  is  the  contact 
resistance  between  the  sample  and  the  holder,  AT  is  the  tempera¬ 
ture  drop  across  the  sample,  and  qs  is  the  heat  transferred  through 
the  sample  which  is  given  by  Eq.  (4).  The  thermal  contact  resistance 
is  calculated  from  the  total  resistance  for  a  single  GDL  as  follows: 

^total  =  ^sample  3"  -sample  (6) 

Since  it  is  not  feasible  to  estimate  the  contact  resistance  be¬ 
tween  the  sample  and  the  holder  using  a  single  sample  with  a  given 
thickness,  we  need  either  to  change  the  thickness  of  the  sample 
used  or  use  a  stack  of  samples.  However,  it  is  significantly  difficult 
to  vary  the  thickness  of  the  sample  as  typically  the  GDL  is  not 
available  in  different  thicknesses.  Alternatively,  the  thickness  is 
varied  by  making  stacks  of  different  numbers  of  individual  GDL 
samples  and  assuming  that  the  contact  resistance  between  the 
samples  are  the  same  [11  ].  However,  this  introduces  a  new  variable, 
which  is  the  contact  resistance  between  individual  samples.  So,  the 
total  contact  resistance  for  a  stack  of  GDL  samples  is  calculates  as 
follows: 


^total,n  —  ^sample  "T  (tl  1  )^sample-sample  T-  ^^H-sample 


In  the  current  study,  a  stack  of  five  GDL  samples  is  used  to  vary 
the  thickness  of  the  tested  specimens.  This  means  that  there  will  be 
five  equations  with  three  unknowns  in  order  to  ensure  that  the 
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Fig.  2.  The  LabVIEW  programme  interface. 
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Fig.  3.  A  typical  steady  state  temperature  profile  through  the  fixture. 


assumption  that  the  contact  resistance  between  the  samples  are 
the  same. 

The  thermal  resistance  offered  by  the  sample  is  given  by: 

^total,l  =  ^sample  T-  2Rj-[_sample  (8) 

The  thermal  resistance  offered  by  the  sample  is  given  by: 


Ls 

ks^s 


where  Ls  is  the  length  of  the  sample,  ks  is  the  thermal  conductivity, 
and  y4s  is  the  heat  transfer  area  [12  .  For  a  given  compression,  the 
contact  resistance  between  the  GDL  sample  and  the  holder  surface 
is  assumed  to  be  insensitive  to  the  number  of  the  samples  involved 
in  the  stack.  Then,  the  thermal  conductivity  of  the  sample  ks  is 
calculated  from  the  slope  of  the  curve  of  the  total  thermal  resis¬ 
tance  as  a  function  of  GDL  thickness: 


1  1 
slope  As 


(10) 


2.5.  Validation  of  the  measurement  technique 


The  uncertainty  in  the  through-plane  thermal  conductivity 
value  is  calculated  and  reported  for  all  the  measurements  and  it 
was  not  more  than  8%  in  this  study. 

3.  Results  and  discussions 

3.1.  Effect  of  the  temperature  on  the  through-plane  thermal 
conductivity  of  the  membrane 

As  seen  in  Fig.  4a,  the  through-plane  thermal  conductivity  of  the 
membrane  is  in  good  agreement  with  the  in-plane  thermal  con¬ 
ductivity  of  the  membrane,  which  has  been  previously  measured  by 
using  the  parallel  thermal  conductance  technique  which  is  based  on 
the  thermal  resistance  measurements  from  the  parallel  resistance 
circuit  depends  on  the  output  power  and  the  temperature  drop 
through  the  sample  [17].  The  maximum  through-plane  thermal 
conductivity  of  the  membrane  was  0.193  ±  0.02  W  m1  K-1  when  the 
temperature  was  35  °C,  and  as  the  temperature  increases  then  the 
thermal  conductivity  of  the  membrane  decreases.  For  example,  the 
minimum  through-plane  thermal  conductivity  of  the  membrane 
was  0.132  ±  0.02  W  m1  I<-1  when  the  temperature  was  80  °C. 

This  decrease  in  the  thermal  conductivity  of  the  membrane  is 
due  to  the  increase  in  the  thermal  resistance  of  the  membrane  with 
temperature,  as  shown  in  Fig.  4b.  Basically,  the  thermal 
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The  experimental  technique  was  validated  by  using  a  standard 
thin  film  material.  The  thermal  conductivity  of  a  stack  of 
27  ±  2.5  pm  aluminium  foils  was  measured  to  be  about 
236.3  ±  4.2  W  m-1  IC1,  which  was  found  to  be  in  good  agreement 
with  the  reported  value  for  the  thermal  conductivity  of  aluminum 
foils,  which  is  about  235  W  m-1  K-1  at  room  temperature  15]. 

2.6.  Uncertainty  analysis 

The  uncertainty  in  the  thermal  conductivity  values  can  be 
estimated  based  on  the  combined  uncertainties  in  the  dimensions 
of  the  sample,  and  the  temperature  drop  across  the  sample  16]  as 
follows: 


where  Tc  and  Th  are  the  temperature  of  the  cold  and  hot  plates, 
respectively.  L  is  the  length  of  the  sample  and  A  is  the  sample  cross 
sectional  area. 


Fig.  4.  Measured  thermal  conductivity  and  thermal  resistance  of  the  membrane  as  a 
function  of  the  temperature  along  with  the  experimental  error  bars. 
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Fig.  5.  Measured  thermal  resistances  of  the  GDLs  as  a  function  of  the  temperature  for  (a)  SGL  10AA,  (b)  SGL  10BA,  (c)  SGL  10CA,  (d)  SGL  10DA,  and  (e)  SGL  10EA. 


conductivity  of  the  membrane  decreases  with  increasing  temper¬ 
ature  due  to  the  decrease  in  the  phonon  mean  free  path  as  a  result 
of  increasing  the  number  of  phonon  at  the  high  temperature  for  the 
polymers  as  Nation  material  [18  . 

This  relation  is  explained  by  using  this  equation  19]: 


^mem  —  ^  ^ph^ph^ph  0 2) 

where  /<mem  is  the  thermal  conductivity  of  the  membrane,  Jph  is  the 
phonon  mean  free  path,  cph  is  the  phonon  heat  capacity  and  vph  the 
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phonon  velocity.  In  the  high  temperature  cph  and  vph  are  almost 
constants,  while  is  JPh  decreases  and  consequently  the  thermal 
conductivity  of  the  membrane  decreases  [20,21],  which  is  consid¬ 
ered  as  an  insulating  material. 

3.2.  Effect  of  the  temperature  on  the  through-plane  thermal 
conductivity  of  the  GDL 

Due  to  the  presence  of  the  PTFE  binder  and  the  polymeric  resin, 
whose  thermal  conductivities  decrease  with  increasing  tempera¬ 
ture,  the  total  thermal  resistance  of  the  GDL  was  found  to  increase 
with  increasing  temperature  and  subsequently  the  thermal  con¬ 
ductivity  decreases  as  the  temperature  increases  [22].  This  effect  is 
clearly  shown  in  Fig.  5. 

3.3.  Effect  of  the  PTFE  loading  on  the  through-plane  thermal 
conductivity  of  the  GDLs 

The  through-plane  thermal  conductivity  of  GDL  is  significantly 
lower  than  its  in-plane  thermal  conductivity  17].  For  example,  at 
35  °C,  the  in-plane  thermal  conductivity  of  10AA  was  found  to  be 
about  12.5  W  m-1  K_1  [17]  whereas  the  through-plane  thermal 
conductivity  of  the  same  material  at  the  same  temperature  was 
measured  as  0.55  W  m”1  K-1  at  2  bar.  It  should  be  noted  that  the 
compression  has  a  positive  effect  on  the  thermal  conduction,  see 
Section  3.4. 

Since  the  heat  transfer  is  from  fibre  to  fibre  in  the  through-plane 
direction,  the  addition  of  PTFE  decreases  the  through-plane  thermal 
conductivity  of  the  gas  diffusion  layer  because  of  the  low  thermal 
conductivity  of  the  PTFE,  which  is  about  0.25  W  m1  K"1  [23-25]. 

As  illustrated  in  Fig.  6,  overall,  the  thermal  conductivity  of  the 
GDL  significantly  decreases  after  PTFE-treatment.  The  GDL,  without 
PTFE,  has  the  highest  through-plane  thermal  conductivity  and  its 
thermal  conductivity  is  about  50%  higher  than  that  of  the  PTFE- 
treated  GDLs.  However,  the  real  relation  between  adding  PTFE 
and  the  through-plane  thermal  conductivity  of  the  GDL  is  not  clear 
as  there  was  no  significant  difference  in  the  through-plane  thermal 
conductivity  of  the  GDL  after  adding  5, 10,  20,  and  30%  PTFE  to  the 
GDL.  This  point  certainly  requires  further  investigation. 

3.4.  Effect  of  compression  pressure  on  the  through-plane  thermal 
conductivity  of  the  GDL 

In  PEM  fuel  cells,  the  GDL  is  deforming  under  the  ribs  and  the 
thickness  of  the  GDL  under  the  ribs  is  less  than  its  thickness  under 
the  channel.  The  GDL  thickness  is  greatly  affected  by  the  external 
compression,  which  affects  the  transport  properties,  the  contact 
resistance  and  the  thermal  and  electrical  conductivities  of  the  GDL 
and  these  subsequently  affect  the  temperature  and  the  water 
management  in  PEM  fuel  cells  [23  .  Therefore,  the  effect  of  the 
compression  on  the  thermal  conductivity  of  the  GDL  is  needed  to  be 
taken  into  account.  The  effect  of  compression  pressure  on  the 
reduction  in  the  thickness  of  the  GDL  is  shown  in  Fig.  7. 

The  reduction  in  the  thickness  is  initially  very  high  when 
applying  the  compression  pressure.  Then,  the  variation  in  thickness 
is  almost  negligible  after  increasing  the  pressure  which  indicates 
that  the  deformation  of  the  GDL  has  ‘saturated’.  Note  that  there  is  a 
hysteresis  effect  in  the  compression  curves  which  signals  that  there 
has  been  a  permanent  deformation  in  the  compressed  GDL  sample 
[26].  In  Fig.  8a,  the  effect  of  the  applied  load  on  the  thermal  resis¬ 
tance  has  been  investigated.  It  is  clear  that  the  thermal  resistance 
values  of  the  treated  GDLs  are  higher  than  that  of  the  untreated 
GDL.  This  is  due  to  the  increase  in  the  contact  resistance  between 
the  fibres  of  the  GDL  after  adding  the  PTFE  [27]. 


Fig.  6.  Measured  thermal  conductivity  of  the  GDLs  as  a  function  of  the  temperature 
and  PTFE  loading  along  with  the  experimental  error  bars. 


The  through-plane  thermal  conductivity  of  the  GDL  increases 
when  the  compression  pressure  increases,  as  shown  in  Fig.  8b.  This 
is  as  a  result  of  the  larger  contact  area  and  consequently  this  pro¬ 
vides  more  heat  paths  and  thus  the  higher  thermal  conductivity 
[28].  It  should  be  noted  that  these  results  are  those  of  increasing 
compressive  pressure. 


3.5.  Effect  of  the  MPL  coating  on  the  through-plane  thermal 
conductivity  of  the  GDL 

In  order  to  enhance  the  water  and  thermal  management  in  the 
PEM  fuel  cells,  a  micro  porous  layer  (MPL)  is  added  to  the  gas 
diffusion  layer.  This  layer  consists  of  a  carbon  black  powder  and 
hydrophobic  agent  [8]. 

The  thermal  contact  resistance  of  the  GDL  samples  decreases 
when  the  MPL  is  added  to  the  GDL.  This  is  not  surprising  as  one  of 
the  main  aims  of  adding  MPL  to  the  GDL  in  the  fuel  cell  is  to 
improve  the  contact  between  the  GDL  and  the  catalyst  layer. 

Furthermore,  it  is  clear  from  Fig.  9  that  the  untreated  10BA  has  a 
lower  thermal  conductivity  than  that  of  10BC  and  10BE,  which  are 
MPL  coated  GDLs.  The  main  reason  for  this  is  that  the  MPL  is  rich  in 
carbon  powder  which  is  considered  as  a  high  thermal  conductivity 
material  compared  to  the  thermal  conductivity  of  the  untreated  GDL. 


3.6.  Effect  of  the  temperature  on  the  through-plane  thermal 
conductivity  of  the  catalyst  layer 


The  through-plane  thermal  conductivity  of  the  catalyst  layer 
was  investigated  by  manually  spraying  the  Pt  ink  onto  some  10BA 
GDL  samples  and  then  measure  the  through-plane  thermal  con¬ 
ductivity  of  the  samples  before  and  after  adding  the  catalyst  layer. 
Since  the  layers  are  lumped  in  series,  the  series  model  [29,30] 
illustrated  in  Fig.  10  was  used  to  calculate  the  thermal  conductiv¬ 
ity  of  the  catalyst,  as  follows: 


1 

Vcl  I  VGDL 

l<  cl  1<G  DL 


(13) 


where  /<totai  is  the  thermal  conductivity  of  the  catalysed  GDL.  /<Gdl 
and  /<d  are  the  thermal  conductivities  of  the  GDL  and  catalyst  layer, 
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Fig.  7.  Hysteresis  in  the  thickness  of  the  GDL  under  the  external  compression  for  (a)  SGL  10AA,  (b)  SGL  10BA,  (c)  SGL  10CA,  (d)  SGL  10DA,  and  (e)  SGL  10EA. 


respectively,  j'gdl  and  vc\  are  the  volume  fractions  of  the  GDL  and 
catalyst  layer,  respectively. 

As  shown  in  Fig.  11,  the  through-plane  thermal  conductivity  of 
the  catalyst  layer,  with  0.4  mg  cm-2  platinum  loading,  was  found  to 


be  comparable  with  the  in-plane  thermal  conductivity  of  the 
catalyst  layer  which  has  been  previously  measured  and  reported  in 
Ref.  [17  .  This  indicates  that  the  catalyst  layer  is  a  homogenous 
material.  Furthermore,  the  through-plane  thermal  conductivity  of 
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Fig.  10.  Calculating  the  thermal  conductivity  of  the  catalyst  using  the  series  model. 
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Fig.  11.  Measured  thermal  conductivity  of  the  catalyst  layers  as  a  function  of  the 
temperature  along  with  the  experimental  error  bars  (the  platinum  loading  for  this  case 
is  0.4  mg  cm-2). 


Fig.  8.  Measured  thermal  conductivity  and  thermal  resistance  of  the  GDLs  at  different 
compression  loads. 


the  catalyst  layer  was  almost  insensitive  to  the  temperature.  The 
maximum  through-plane  thermal  conductivity  of  the  catalyst  layer 
is  0.345  ±  0.02  W  m-1  I<-1  when  the  temperature  is  35  °C,  and  the 
minimum  through-plane  thermal  conductivity  of  this  catalyst  layer 
is  0.334  ±  0.02  W  m-1  I<-1  when  the  temperature  is  80  °C. 

In  order  to  investigate  the  effects  of  the  Pt  loading  on  the 
through-plane  thermal  resistance  of  the  catalyst  layer,  the  catalyst 


layers  with  three  different  Pt  loadings  (0.2,  0.4,  0.6)  mg  cm-2  were 
prepared  and  the  through-plane  thermal  resistance  of  them  were 
compared,  see  Fig.  12. 

The  through-plane  thermal  resistance  of  the  catalyst  layer  was 
found  to  decrease  with  Pt  loading.  This  is  mainly  because  of  the 
thickness  of  the  catalyst  layer.  As  the  platinum  loading  increases, 
the  thickness  of  the  catalyst  layer  increases.  This  leads  to  a  better 
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Fig.  9.  Measured  thermal  conductivities  for  (a)  SGL  10BA,  (b)  SGL  10BC,  and  (c)  SGL 
10BE  along  with  the  experimental  error  bars. 


Fig.  12.  Measured  thermal  resistances  of  the  catalyst  layers  as  a  function  of  the  Pt 
loading  along  with  the  experimental  error  bars. 
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contact  between  the  catalyst  layer  and  the  holder  of  the  sample. 
Therefore,  the  reduction  in  the  total  resistance  as  a  result  of 
increasing  the  platinum  loading  is  more  to  do  with  the  contact 
resistance  and  not  with  the  thermal  conductivity  which  must 
ideally  remain  constant. 

4.  Conclusions 

In  this  study,  an  experimental  technique  has  been  developed  to 
determine  the  through-plane  thermal  conductivity  and  the  contact 
resistance  of  the  components  in  the  MEA  and  GDL.  The  main  con¬ 
clusions  of  these  measurements  are  as  follows: 

•  The  through-plane  thermal  conductivity  of  the  membrane  was 
found  to  be  comparable  with  the  in-plane  thermal  conductivity 
and  it  decreases  with  increasing  the  temperature.  The  through- 
plane  thermal  conductivity  of  the  membrane  was  found  to  be 
0.193  ±  0.018  W  m-1  K-1  and  0.132  ±  0.02  W  irr1  K"1  when  the 
temperature  was  35  °C  and  80  °C,  respectively. 

•  The  through-plane  thermal  conductivity  of  the  GDL  was  found 
to  be  significantly  lower  than  its  in-plane  thermal  conductivity. 
The  through-plane  thermal  conductivity  of  the  GDL  decreases 
with  increasing  the  temperature  and  increases  with  increasing 
compression  pressure. 

•  The  through-plane  thermal  conductivity  of  the  catalyst  layer  is 
comparable  to  its  in-plane  thermal  conductivity  and  it  was 
found  to  be  almost  insensitive  to  the  temperature  and  the 
thermal  resistance  of  the  catalyst  layer  decreases  slowly  with  Pt 
loading  due  to  the  increase  in  the  thickness  of  the  catalyst  layer 
and  the  contact  resistance.  The  through-plane  thermal  resis¬ 
tance  of  the  catalyst  layer,  when  the  temperature  was  35  °C,  was 
found  to  be  1.985  ±  0.18  K  W” 1  and  1.553  ±  0.13  K  W”1  when  the 
Pt  loading  was  0.2  and  0.6  mg  cm'2,  respectively. 

The  through-plane  thermal  conductivity  of  the  GDL  and  the 
MEA  components  provides  some  of  the  necessary  parameters 
which  will  enhance  the  predictions  of  the  profiles  of  temperature 
and  water  saturation  in  PEM  fuel  cells.  As  a  future  work,  it  will  be  of 
great  interest  to  measure  the  thermal  conductivity  of  the  compo¬ 
nents  of  the  MEA  beyond  the  temperature  of  120  °C,  which  is 
normally  the  operating  temperature  for  high-temperature  PEM  fuel 
cells. 
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Nomenclature 

A  cross-sectional  area  of  the  sample  (m2) 

cph  the  phonon  heat  capacity  (J  g-1  K_1) 
ks  thermal  conductivity  of  the  sample  (W  m-1  I<-1) 

km  thermal  conductivity  of  the  standard  material 

(W  m  1  K"1) 

/{mem  thermal  conductivity  of  the  membrane  (W  m_1  K-1) 
kc\  thermal  conductivity  of  catalyst  layer  (W  m-1  K-1) 

/<gdl  thermal  conductivity  of  GDL  (W  m-1  I<-1) 

/Ph  the  phonon  mean  free  path(m) 

Ls  length  of  the  sample  (m) 


qiow  heat  transfer  through  the  lower  plate  (W) 

qs  heat  transfer  through  the  sample  (W) 

qup  heat  transfer  through  the  upper  plate  (W) 

^H-Sampie  thermal  contact  resistance  between  two  layers  of  the 
sample  (m2  KW^1) 

^sample-sample  thermal  contact  resistance  between  two  layers  of  the 
sample  (m2  KW-1) 

Rs  sample  thermal  resistance  (K  W_1) 

i^iotai  total  thermal  resistance  (I<  W-1) 
vph  the  phonon  velocity(m  s-1) 

T\  temperatures  measured  by  the  first  thermocouple  (°C) 

T3  temperatures  measured  by  the  third  thermocouple  (°C) 

T4  temperatures  measured  by  the  fourth  thermocouple  (°C) 

Tq  temperatures  measured  by  the  sixth  thermocouple  (°C) 

7c  cold  plate  temperature  (°C) 

7h  hot  plate  temperature  (°C) 

A T  temperature  drop  (K) 

vc\  volume  fraction  of  catalyst 

J'gdl  volume  fraction  of  GDL 

Abbreviations 

CL  catalyst  layer 

GDL  gas  diffusion  layer 

MEA  membrane  electrode  assembly 

MPL  micro  porous  layer 

PEM  proton  exchange  membrane 

PTFE  polytetrafluoroethylene 
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